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ABSTRACT 

Determining the temperature distribution of coronal plasmas can provide stringent constraints on 
coronal heating. Current observations with the Extreme ultraviolet Imaging Spectrograph onboard 
Hinode and the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory provide 
diagnostics of the emission measure distribution (EMD) of the coronal plasma. 

Here we test the reliability of temperature diagnostics using 3D radiative MHD simulations. We 
produce synthetic observables from the models, and apply the Monte Carlo Markov chain EMD 
diagnostic. By comparing the derived EMDs with the "true" distributions from the model we assess 
the limitations of the diagnostics, as a function of the plasma parameters and of the signal-to-noise of 
the data. 

We find that EMDs derived from EIS synthetic data reproduce some general characteristics of the 
true distributions, but usually show differences from the true EMDs that are much larger than the 
estimated uncertainties suggest, especially when structures with significantly different density overlap 
along the line-of-sight. When using AIA synthetic data the derived EMDs reproduce the true EMDs 
much less accurately, especially for broad EMDs. The differences between the two instruments are 
due to the: (1) smaller number of constraints provided by AIA data, (2) broad temperature response 
function of the AIA channels which provide looser constraints to the temperature distribution. 

Our results suggest that EMDs derived from current observatories may often show significant dis- 
crepancies from the true EMDs, rendering their interpretation fraught with uncertainty. These in- 
herent limitations to the method should be carefully considered when using these distributions to 
constrain coronal heating. 

Subject headings: X-rays, Sun, EUV, spectroscopy; Sun: corona 



1. INTRODUCTION 

The heating mechanism that is responsible for the mil- 
lion degree solar corona remains unknown, though sev- 
eral candidates exist. It is one of the most important 
open issues in astrophysics. Constraining the properties 
of this heating mechanism is a difficult task, but usually 
performed through spectral and imaging observations of 
the solar corona, which provide diagnostics of the plasma 
temperature distribution. The latter has important im- 
plications for the energy balance of the corona (see e.g., 
Klimchuk 2006; Reale 2010 and references therein). 

The thermal distribution of the plasma - or emission 
measure distribution, EMD (in section 2 we define the 
relationship between EMD and the differential emission 
measure, DEM, which is also often used to describe the 
plasma thermal distribution) - is crucial to test heating 
models. For instance, the EMD in coronal loops, strongly 
depends on the spatial and temporal properties of the 
energy release (e.g., Klimchuk & Cargill 2001; Cargill & 
Klimchuk 2004; Testa et al. 2005). The presence of a high 
temperature (T > 5MK) component, e.g., in the EMD 
of an active region, is a good tracer of the properties of 
the heating (e.g., Patsourakos & Klimchuk 2006), and 
it has recently been addressed by several studies (e.g., 
Patsourakos & Klimchuk 2006; Reale et al. 2009b, a; Ko 
et al. 2009; Schmelz et al. 2009a; Testa et al. 2011; Testa 
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& Reale 2012). Also, the thermal distribution in the 
cross-field direction can help discern between a mono- 
lithic, single strand, loop (isothermal at a given location 
along the loop; e.g., Aschwanden et al. 2000; Aschwan- 
den & Nightingale 2005; Aschwanden & Boerner 2011; 
Del Zanna 2003; Landi et al. 2002, 2006; Landi & Feld- 
man 2008), and a loop structure composed of several 
strands (multi-thermal plasma; e.g., Schmelz et al. 2001, 
2005; Warren & Brooks 2009; Brooks et al. 2009). 

Ample efforts have been devoted to the accurate de- 
termination of the thermal structuring of coronal plasma 
to derive robust observational constraints on the coronal 
heating mechanism(s). The plasma temperature distri- 
bution of the quiet corona and of active regions has been 
investigated through imaging data and spectroscopic ob- 
servations (e.g., Brosius et al. 1996; Landi & Landini 
1998; Aschwanden et al. 2000; Testa et al. 2002; Del 
Zanna & Mason 2003; Reale et al. 2007; Landi et al. 2009; 
Shestov et al. 2010; Sylwester et al. 2010). Several re- 
cent studies have focused on EUV spectra obtained with 
the Hinode Extreme Ultraviolet Imaging Spectrometer 
(EIS; Culhane et al. 2007) which provides good temper- 
ature diagnostic capability, together with higher spatial 
resolution and temporal cadence than previously avail- 
able (e.g., Watanabe et al. 2007; Warren et al. 2008; 
Patsourakos & Khmchuk 2009; Brooks et al. 2009; War- 
ren & Brooks 2009; Testa et al. 2011; Tripathi et al. 
2011). Imaging observations obtained by the SDO At- 
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Fig. 1. — Three-dimensional snapshot of one of the 3D rMHD models (model H): top view (left panel) and side view (right panel). Some 
selected magnetic field lines are shown in red to give an indication of the magnetic field topology in the corona. The isosurfaces of Fe xil 
and Caxvil emission are shown in green and blue respectively. We also show at the photospheric layer (z=0 Mm) the strength of the 
vertical magnetic field in blue-red color scale. 



mospheric Imaging Assembly (AIA; Lemen et al. 2012) 
in narrow EUV bands have also been used for studying 
the plasma temperature distribution, especially in active 
region loops (e.g., Schmelz et al. 2010; Aschwanden & 
Boerner 2011; Aschwanden et al. 2011; Schmelz et al. 
2011), also in conjunction with EIS data (e.g.. Brooks 
et al. 2011; Warren et al. 2011; Testa & Reale 2012). 

The optically thin nature of the coronal EUV and 
X-ray emission implies that the emission observed in 
the resolution element (instrument pixel) is generally 
produced by several independent structures along the 
LOS. Most plasma diagnostics rely on some homogene- 
ity assumptions (e.g., constant plasma density along the 
LOS) whereas these overlapping structures can in princi- 
ple have very different plasma parameters. This poses 
significant challenges for interpreting the meaning of 
the derived plasma parameters (which are by necessity 
weighted averages of the distributions along the LOS) in 
terms of the actual physical conditions of the plasma. 
Further assumptions typically made for specific diagnos- 
tics (e.g., on the functional form and smoothness of the 
EMD) also have a potentially significant impact on the 
accuracy of the diagnostics. Several efforts have been 
carried out to test the accuracy of the plasma tempera- 
ture diagnostics. For instance, the notorious challenges 
in determining the emission measure distribution and its 
confidence limits have been addressed in several studies 
(e.g., Craig & Brown 1976; Judge et al. 1997; Mcintosh 
2000; Judge 2010; Landi & Klimchuk 2010; Landi et al. 
2012; Hannah & Kontar 2012). The main limitation of 
these previous efforts lies in the necessarily simplified test 
cases adopted. 

In this paper we use advanced 3D radiative MHD sim- 
ulations of the solar atmosphere with the Bifrost code 
(Gudiksen et al. 2011) to carry out detailed tests of 
coronal plasma temperature diagnostics. We focus on 
the plasma temperature diagnostics using current spec- 
tral (Hinode/ElS) and imaging {SDO/AIA) data. We 
synthesize EIS and AIA data from 3D radiative MHD 
(rMHD) simulations and analyze them like real data, and 
use the comparison of the thermal distributions inferred 
from the synthetic data with the "true" distributions, 
which in the case of the simulations are known, to care- 
fully assess the accuracy of the diagnostics. These 3D 
simulations provide us with the opportunity to improve 



upon previous work by exploring more realistic configu- 
rations, with significant superposition of different struc- 
tures along the LOS, allowing a statistical approach to 
determine the accuracy and limitations of the plasma 
diagnostics, for a variety of spatial and thermal structur- 
ing of the plasma. The three dimensional nature of the 
simulations also allows us to explore a variety of realis- 
tic viewing angles, reproducing typical distributions of 
structuring ranging from on disk to limb observations. 

In Section 2 we describe the analysis methods, we in- 
clude a short description of the Bifrost code, and discuss 
the characteristics of the 3D rMHD models and of the 
corresponding synthetic observables used in this work. 
The analysis of the synthetic data and the results of the 
determination of the plasma temperature distribution are 
presented and discussed in Section 3. We summarize our 
findings and draw our conclusions in Section 4. 

2. ANALYSIS METHODS 

In this work we investigate how accurately the thermal 
distribution of the plasma can be inferred from coronal 
observations that are currently available. The observed 
intensities of a set of spectral emission lines constrain 
the plasma temperature distribution, as they depend on 
the abundance Az of the emitting element, the plasma 
emissivity of the spectral feature GxiT^Ue) as a func- 
tion of temperature T and electron density ng, and the 
differential emission measure distribution DEM(T)\ 

Ix=Az Gx(T, n,)DEM(T) dT (1) 

where DEM{T) = nldV/dT [cm-^R-^]. Analogously, 
in the case of imaging observations in broad/narrow pass- 
bands the observed intensity in a channel will depend on 
the temperature response function, i^chan(^)* 

/chan = ^ RcY..n{T)DEM{T) dT. (2) 

Throughout the paper, instead of the DEM(T) we will 
discuss the emission measure distribution EM(T) which 
is obtained by integrating the differential emission mea- 
sure distribution DEM{T) in each temperature bin; here 
we use a temperature grid with constant binning in log- 
arithmic scale (AlogT = 0.05). 
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Fig. 2. — Histogram of the electron density (left panel) and temperature (right panel) for the two used snapshots (for the voxels with 
logT > 5.3). The simulation "C" modeling the smaller region with emerging flux (solid line) is characterized by lower coronal temperature 
than the larger simulation "H" (dashed lines), and it has a broader distribution of densities (where the tail at logne ^ 10 corresponds to 
dense plasma in the emerging flux region). 



Several methods have been developed to reconstruct 
emission measure distributions from a set of observed 
intensities in lines, or passbands in the case of imag- 
ing observations (see e.g., review by Phillips et al. 
2008, and the discussion and references in Hannah 
& Kontar 2012). Here we test the Monte Carlo 
Markov chain (MCMC, hereafter) forward modeling 
method (Kashyap & Drake 1998), which is widely 
used and considered to provide robust results (see e.g., 
Landi et al. 2012; Hannah & Kontar 2012; see also 
http://www.lmsal.com/~boerner/demtest/ for a recent 
comparative analysis of results from different methods, 
applied to AIA). With respect to several other methods, 
the MCMC method has the advantages of not impos- 
ing a pre-determined functional form for the solution, 
and, most importantly, of estimating the uncertainties 
associated with the resulting emission measure distribu- 
tion (see e.g., Kashyap & Drake 1998; Testa et al. 2011, 
for additional details). We use the Package for Inter- 
active Analysis of Line Emission (PINTofALE, Kashyap 
& Drake 2000) which is available as part of Solar Soft. 
Though no pre-determined functional form is imposed, 
the MCMC method does apply some smoothness criteria 
which are locally variable and based on the properties of 
the temperature responses/emissivities for the used data, 
instead of being arbitrarily determined a priori. 

For a given simulation snapshot, which provides the 
electron density Uq and temperature values for each grid 
point (voxel) in a three dimensional box (see section 2.1 
for a description of the characteristics of the simulations 
we analyzed), and a selected LOS we proceed as follows: 

• Using the Uq and T values, and the emissivities 
from CHIANTI (Dere et al. 1997, 2009) we synthe- 
size intensities of a set of EIS lines and in the 6 
AIA coronal channels, using the optically thin ap- 
proximation and statistical ionization equilibrium; 

• by integrating through the box along the LOS, and 
degrading the spatial resolution to the instrument 
resolution, we obtain synthetic coronal images in 
the different lines/passbands; 

• we consider two cases: with or without photon 
(Poisson) noise; i.e., in the latter case we random- 
ize the intensities according to the photon count- 



TABLE 1 

EIS LINES SYNTHESIZED FROM 3D MODELS 
AND USED FOR RECONSTRUCTING THE 
EMISSION MEASURE DISTRIBUTION. 



A[A] ^ 


Ion 


log(Tn.ax[K]) 


Notes ^ 


268.991 


Mgvi 


5.65 


sb 


185.213 


Fe VIII 


5.70 




278.404 


Mgvii 


5.80 




275.361 


Si VII 


5.80 




188.497 


Fe IX 


5.90 




184.537 


Fex 


6.05 




188.216 


Fexi 


6.15 




195.119 


FeXH 


6.20 


sb 


274.204 


Fexiv 


6.30 




284.163 


Fexv 


6.35 




262.976 


Fexvi 


6.45 




208.604 


Caxvi 


6.70 




192.853 


Caxvii 


6.75 


bl 


254.347 


Fe XVII 


6.75 





^ The wavelengths are from CHIANTI (in 
case of self-blend we list the wavelength of 
the strongest line). 

^ The label "sb" indicates that the spectral 
feature is a self-blend of lines from the same 
ion, all included in synthesizing the data. 
The label "bl" for the Caxvii line indicates 
that at the EIS spectral resolution this line 
is blended with Fe xi and O V lines (Ko et al. 
2009). For this paper however we only syn- 
thetize the Caxvii intensity, i.e., we do not 
calculate the blending lines. 

ing statistics. The noise level and the uncertain- 
ties associated with the intensities are calculated 
by assuming signal-to-noise ratios typical of actual 
observations; 

• we use the intensities and uncertainties, in each 
pixel, as input for the MCMC routine to calculate 
the EM{T) solution, pixel by pixel; 

• we analyze the results, using several parameters to 
assess the ability of the method to reproduce the 
input intensities and the "true" EM(T)^ which in 
these test cases are known. 

In the following subsection, section 2.1, we describe in 
detail the choices and assumptions made, and the char- 



4 



Testa et al. 




acteristics of the selected simulations and of the resulting 
synthetic data. 

2.1. 3D Models and Synthetic Observables 

The model considered spans from the upper layer of 
the convection zone up to the low corona, and self- 
consistently produces a chromosphere and hot corona, 
through Joule dissipation of electrical currents. In order 
to investigate the robustness of the EM(T) reconstruc- 
tion method for a wide range of plasma conditions, we 
selected snapshots from two simulations that are char- 
acterized by significantly different parameters (strength 
and spatial distribution of the seed magnetic field, and 
dimensions of the box) , which leads to significantly differ- 
ent distributions of temperature and density throughout 
the box. 

We analyze one of the snapshots from the simulation 
that was previously used to investigate the relative con- 
tribution of different lines to the emission in the AIA 
passbands (Martmez-Sykora et al. 2011). This simula- 
tion covers a volume {x^y^z^ where z is the vertical di- 
rection) of 16 X 8 X 16 Mm^ (512 x 256 x 365 grid points; 
with Ax = Ay ~ 32km, and Az non-uniform and smaller 
where gradients are large with values Az ~ 28km up to 
a height of 4Mm and increasing to « 150km at the top of 
the computational box), and is highly dynamic yielding 
coronal temperatures that increase with time. We used 
an intermediate snapshot {t = 1200s), where the distri- 
bution of temperature peaks around log{T[K]) ~ 6.25, 



and logTmax ^ 6.4. As we discussed in Martmez-Sykora 
et al. (2011), the plasma conditions of this simulation 
represent a good comparison for quiet Sun/coronal hole 
conditions, with some small emerging flux regions, e.g., 
small bright points. In the following we use the label C 
to refer to this cooler and smaller snapshot. 

The other snapshot we consider (H, hereafter), is from 
a simulation modeling a larger region (24x24x 16Mm^; 
768 X 768 X 768 grid points; Ax = Ay 31km, Az ^ 
14km up to a height of 4Mm and increasing to ~ 80km 
at the top of the computational box; see Figure 1), with 
a magnetic field configuration consisting of two small re- 
gions of opposite polarities similar to a small active re- 
gion (Carlsson et al., in preparation). The distribution 
of coronal temperatures is also similar to typical active 
region values, with a peak around log(T[i^r]) ~ 6.6, and 
logTmax ^ 6.7. 

Figure 2 shows the distributions of electron density 
and temperature throughout the simulation box, for the 
two selected snapshots, including only plasma with logT 
higher than 5.3. These distributions show that simula- 
tion "C" , modeling the smaller region with emerging flux, 
is characterized by lower coronal temperature than the 
larger simulation "H" , and it has a broader distribution 
of densities. 

Intensities of a selection of EIS lines and AIA pass- 
bands were computed for each voxel of the 3D simula- 
tion box Vxyz^ using the plasma electron density and 
temperature values and using the density and temper- 
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Fig. 4. — Synthetic top- view (top row) and side-view (bottom row) Fexil (195.1 A) and Caxvil (192. SA) images, at the intrinsic resolution 
of the simulation, from the larger and hotter snapshot (H) of the numerical simulation containing a plage-like region, representative of the 
core of a small active region. 



ature dependent contribution functions G(T, ng) from 
CHIANTI (Dere et al. 1997, 2009), assuming coronal 
abundances (Feldman 1992), and the ionization fractions 
of chianti.ioneq. For AIA we compute synthetic intensi- 
ties in the coronal passbands (94A, 131A, 171A, 193A, 
211 A, 335 A), by calculating the spectra (line by line, tak- 
ing into account the temperature and density dependence 
of the contribution function for each line, and adding the 
continuum emission) and folding them through the effec- 
tive area of each channel (Boerner et al. 2012). In the 
case of EIS we selected a set of lines that provide a good 
coverage of the temperature range logT ~ 5.6 — 6.7, but 
using a relatively limited number of lines, representative 
of a typical EIS observation: we selected 14 lines, which 
are listed in Table 1. 

Images are then derived for each snapshot, and for two 
different LOS, by integrating the emission through the 
simulated box along the LOS: along the horizontal direc- 
tion y for the "side view", and along the vertical direc- 
tion z in the case of the "top view" , analogously to the 
case presented in Martmez-Sykora et al. (2011). Figures 
3 and 4 show examples of such images, at the intrinsic 
resolution of the simulations. In addition, we degrade 
the spatial resolution to the instrumental spatial resolu- 
tion: 0.6 arcsec/pixel for the AIA synthetic data, and 
1 arcsec/pixel for EIS. 

In order to compute the uncertainties of the simulated 
intensities we assume a signal-to-noise ratio (S/N) typi- 
cal of well exposed active region observations (e.g., Reale 
et al. 2011). For the bright channels of AIA (171A, 
193A, 211A) we use the DN pix ^ assuming typical ex- 
posure times in a single image (texp ^ 2s). For the other 
channels (94A, 131A, 335A) which typically have signifi- 
cantly fainter emission we use the S/N calculated assum- 
ing the exposure of 10 summed images (corresponding to 
^exp ^ 30s) . In the case of EIS, we calculate the S/N by 
scaling the images so that the Fexil 195A line, which is 
usually one of the brightest emission lines in EIS obser- 



vations, has 500 counts, as in typical observations, and 
computed the uncertainties from photon counting statis- 
tics. We first considered the case without Poisson noise, 
and then included the Poisson noise by randomizing the 
intensity values assuming the photon counting statistics 
of typical AIA and EIS observations as described above. 

In appendix A we present a comparison of the synthetic 
intensities derived from the simulations, as described 
above, with measured intensities from recent SDO/AIA 
observations. This comparison shows that the simula- 
tions here considered yield emission values of similar or- 
der of magnitude of real observations, and therefore pro- 
vide a sensible test case for the coronal diagnostics. 

3. RESULTS 

Following the procedure described in detail in sec- 
tion 2, we analyzed the synthetic data to diagnose the 
plasma temperature distributions, pixel by pixel, using 
the MCMC reconstruction method. We will first discuss 
the results obtained using the EIS synthetic data, and 
then describe the results obtained with the AIA intensi- 
ties. 

When analyzing real data, the only measure of the 
goodness of the inferred emission measure distribution 
consists in the agreement between the measured inten- 
sities and the emissions predicted using the best fit 
EMD. Therefore, as a first step, we show the maps of 
Xo ( = ^jii^j^pred - Ij^ohs)/crjf /df, whcrc Ij^ohs and (7j 
are the synthetic intensity and associated error for the 
line/channel j, Ij^pred is the intensity predicted by the 
best fit EMD and df are the degrees of freedom^), for 
both snapshots and both LOS, using the EIS synthetic 
intensities (snapshots C and H in Figure 5, and 6, respec- 
tively). We show both the case without Poisson noise 
{left)^ and including Poisson noise (right). In both cases 

^ The degrees of freedom are calculated as de- 
scribed in the PINTofALE documentation http://hea- 
www.harvard.edu/PINTofALE/doc/MCMC_DEM. html#out. 
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Fig. 5. — Maps showing the Xo for the EMD derived using the EIS hnes hsted in Table 1, for snapshot C, for the two LOS, "xy" (top), 
and "xz" (bottom). The right column shows the case where Poisson noise has been included. The crosses mark two pixels with significantly 
different Xo that will be analyzed in more detail later in the paper (see Figure 8). Pixels where the EMD is not calculated - because less 
than 5 lines/channels have non-negligible intensity - are plotted in white color. 



the errors on the intensities, <jj, are calculated as the 
Poisson error on the counts values. For both snapshots 
the Xo maps indicate that the inferred EMD on average 
reproduce the "measured" fluxes adequately. The case 
without noise is characterized by Xo lower than 1 over 
large areas, while areas with large Xo present clear struc- 
turing, which is still present in the case including the 
effect of noise. 

For the top view case, the poorer results in those areas 
appear to be related to the plasma electron density. To 
investigate this effect, we use the diagnostic based on 
the Fexil line ratio 186.88A/195.12A, which is one of 
the most useful density diagnostics accessible with EIS 
spectra (e.g.. Young et al. 2009), and included in many 
EIS observing programs. 

In Figure 7 we show the top view map of electron den- 
sity obtained from the ratio of the synthetic Fe xil inten- 
sities, computed from the two snapshots. By comparing 
these density maps to the Xo plo^s (top row. Figures 5 
and 6), it is clear that the high density regions, such 
as for instance the emerging flux region of snapshot C 
or footpoint ("moss") regions of snapshot H, generally 
correspond to high Xo areas. 

As an example, in Figure 8 we show the electron den- 
sity distribution from the simulations, within two pix- 



els with significantly different Xo (with pixel A showing 
lower Xo than B), in the top view case for each of the 
two snapshots. The selected pixels are marked in Fig- 
ure 5 and 6 (top row panels). We find that the plasma 
in pixel B {xq ^ 2) has a distribution of density along 
the LOS much broader than the plasma in pixel A. This 
leads to a worse determination of the emission measure 
and worse match of the observed fluxes: even if the se- 
lected lines, similar to typical analyses of real data, are 
mostly insensitive to electron density, their small den- 
sity dependency significantly affects the diagnostics of 
emission measure distributions when there is significant 
mixing of plasmas with different electron densities within 
the LOS. The same effect of density mixing is causing the 
large Xo regions at the loop footpoints in the side views 
of both snapshots (bottom panels of Figure 5 and 6). 

For the side view we note that for snapshot H, in 
the higher atmosphere, where the temperature distribu- 
tion in each pixel is essentially isothermal at high tem- 
peratures (see also appendix B), the case including the 
noise tends to have slightly lower Xo values (Figure 6). 
We interpret this effect as due to the extremely nar- 
row (isothermal) distribution: the MCMC method ap- 
plies some physically based, locally variable, smoothness 
criteria based on the properties of the temperature re- 
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Fig. 6. — Xo maps analogous to Figure 5 but for snapshot H. 



sponses/emissivities for the used data (see also discus- 
sion in Kashyap & Drake 1998; Testa et al. 2011), and 
therefore it does not perfectly recover purely isothermal 
distributions; the effect of the noise mimics slight de- 
partures from isothermal distributions and therefore the 
MCMC methods finds better matches to the intensities. 
A similar effect is also observed for small regions of snap- 
shot C (e.g., the high Xo strip at z=19 pix and x from 
pixel 7 to 14). 

We now consider the Xo iiiaps obtained for the analysis 
of the AIA synthetic intensities; Figures 9 and 10 show 
the maps for snapshot C and H respectively (analogous 
to the EIS results of Figures 5, and 6). For snapshot C 
the top view shows a few pixels with higher Xo the 
region corresponding to the dense emerging flux region, 
similarly to the EIS case. For snapshot H there is no clear 
correspondence between the Xo value and the regions 
with high density mixing. This can be explained by look- 
ing at the density sensitivity of the AIA bands, which we 
have already discussed in Martmez-Sykora et al. (2011) 
and summarize in Figure 11. The plots in Figure 11 
show that most channels have limited density sensitivity, 
which becomes negligible at high temperature. There- 



fore we expect less of a density effect for snapshot H 
where high temperature plasma (logT > 6.5; see also 
appendix B) largely contributes to the observed intensi- 
ties. In the side view case, the intensities are reproduced 
quite well in the cases that do not include Poisson noise, 
whereas the cases including noise show large areas of high 
Xo values, especially at large heights. The larger effect 
of the noise for AIA is due to the fact that for the side 
view, the AIA intensities are rather low (see also Fig- 
ure 26) in some channels, especially in 131A for both 
snapshots and 94 A for snapshot C. These channels are 
more sensitive to the cool plasma (though, for snapshot 
H the 94A intensity is largely due to the Fexvill emis- 
sion from the high temperature plasma) which has very 
small or zero contribution higher up in the atmosphere, 
where the plasma is close to isothermal (as we will show 
in detail in section 3.1). The low statistics in few chan- 
nels makes the fit sensitive to the noise. As a result, the 
fits become worse when including the effect of Poisson 
noise. We have investigated the results for significantly 
larger S/N ratio for side view case of snapshot H. We 
re-run the case including noise, but assuming exposure 
times that correspond to 50 images for all AIA channels. 
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Fig. 7. — Maps of electron density derived from synthetic data 
for the two snapshots (C, top, H, bottom), for the top view, using 
the Fexii diagnostics (186A/195A) available with EIS observations 
(see e.g.. Young et al. 2009). 

We find that for very large S/N values the Xo values for 
the intensities tend to get large, because of the very small 
errors associated with the intensities (down to the level 
of ~ 0.1% for the brightest channels). We find that in 
order to get reasonable Xo values the MCMC procedure 
needs to be run with a significantly larger number of 
simulations to reproduce the intensities to that level of 
accuracy. However, as we will discuss later in section 3.2, 
results show that the effect of the noise is not the main 
factor that determines the ability of the MCMC method 
to reconstruct EMDs from AIA data. 

For snapshot H, for both LOS, for the case without 
noise the Xo values are generally worse than the corre- 
sponding C snapshot (Figure 10). We interpret this as 
an effect of the broader temperature distributions (es- 
pecially in the top view) extending to rather high tem- 
peratures {logT[K] > 6.6; see also following Figures 16 
and 17, and appendix B and section 3.1). In this temper- 
ature regime AIA provides less reliable temperature di- 
agnostics than at lower temperatures (log T[i^] ~ 6 — 6.3) 
where several AIA channels have high sensitivity and 
their relative shapes in temperature response provide 
better constraints to the plasma temperature distribu- 
tions. 

We note that the uncertainties we assume here are gen- 
erally significantly smaller than the errors typically as- 
sociated to intensity measurements in real data for the 
EMD reconstruction. Those errors are typically larger 



than the uncertainties associated with photon counting 
statistics because they also take into account uncertain- 
ties in atomic data which are difficult to quantify. In the 
test presented in this paper, there are no uncertainties 
to be associated with the atomic data, since we use the 
same atomic data to synthesize the intensities and to in- 
fer from them the EMD. Nevertheless, the small errors 
we assume can in part explain the sometimes large values 
of Xo we obtained here. 

The results discussed so far provide us with a measure 
of how well the temperature distributions inferred us- 
ing the MCMC method are able to reproduce the "mea- 
sured" intensities. We now discuss in detail the compar- 
ison of the derived EMD with the true emission measure 
distributions. We looked at different parameters to assess 
the robustness and limits of the reconstruction method in 
recovering the plasma temperature distributions: (1) the 
temperature at which the EMD peaks {T {EMD max)) 
within the temperature range log T[i^] = [5.5 — 6.7], (2) 
the EM values integrated in broad temperature bins, and 
(3) the full EMD (i.e., at the fine temperature resolu- 
tion). The temperature of the EMD peak is not nec- 
essarily a good indicator to evaluate the goodness of 
the MCMC solution, especially for broader and more 
structured EMDs. We therefore present the analysis of 
T{EMDmax) in appendix A, summarizing here only the 
main results, while in the main text we will focus on the 
other two comparisons (sections 3.1 and 3.2). 

The results presented in appendix A show that 
the MCMC method performs well in recovering the 
T {EMD max) when the EMD is close to isothermal. For 
broad EMDs the method tends to generally underesti- 
mate T{EMDmax)-, due to several factors, including poor 
constraints beyond the log T[i^r] = [5.5 — 6.7] range, and 
hmited capability of the method to recover sharp fea- 
tures. 

3.1. Emission measure maps 

One of the approaches we use to determine the accu- 
racy of the EMD diagnostics is by looking at the maps 
of emission measure in different temperature ranges. 

In Figure 12 and 13 we show, for top and side view of 
snapshot C respectively, the maps of the true emission 
measure (left column), in six relatively broad tempera- 
ture bins, spanning the temperature range log T[i^] = 
5 — 7. Figures 16 and 17 show the analogous plots for 
snapshot H. Figures 14 and 15 (for snapshot C), and Fig- 
ures 18 and 19 (for snapshot H) also show the EM maps 
derived at the AIA spatial resolution, for direct compar- 
ison with the corresponding maps derived by the analy- 
sis of the synthetic data. These plots show, for the top 
view, how the plasma is characterized by rather broad 
temperature distributions in most locations, especially 
in locations associated with the emerging flux region for 
snapshot C, and with the loop footpoints ("moss") for 
snapshot H (see Figures 3 and 4). The simulated corona 
in snapshot C lacks significant plasma volumes at tem- 
peratures above log T[i^r] = 6.4, while for snapshot H 
the coronal plasma reaches temperature even higher than 
log T[i^] = 6.7 in a hot loop-like feature, which emits 
brightly in e.g., Caxvil (Figure 4). The side view shows 
loop-like features in the lower corona of both snapshots, 
with rather broad temperature distribution. Higher up 
in the corona the plasma is characterized by narrow tem- 
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Fig. 8. — Distribution of electron density, from the simulations, in the voxels contributing to two pixels with different Xo (pixel A with 
Xq < 1; pixel B with Xo ^ 2) in the top view for each snapshot (simulation C in the top row and simulation H in the bottom row) to 
illustrate the effect of the mixing of plasma volumes characterized by significantly different plasma parameters. The selected pixels are 
marked in the plots showing the Xq maps (Figure 5 and 6 for C and H respectively). Histograms of densities in the two pixels are shown 
in the left panel; 2D histograms, as a function of both temperature and densities are also shown (middle: pixel A; right: pixel B). 



perature distributions peaking around log T[K] ~ 6.2 in 
snapshot C, and log T[K] 6.6 in snapshot H. 

The corresponding maps for the EMD derived from the 
analysis of EIS and AIA synthetic data are shown next to 
the true maps for ah the cases, without or including noise 
(in center and right column respectively of Figures 12-15, 
and 16-19). 

For EIS, the general features of the emission measure 
distribution in the central temperature bins (i.e., between 
log T[i^] ~ 5.5 and ~ 6.4 for snapshot C and in the range 
log T[i^] ~ 5.5-6.7 for snapshot H), are recovered quite 
well, also when including the effect of the noise. In the 
top view case, for both snapshots, in the temperature bin 
at the low end of the constrained range (log T[i^r] ~ 5.5 — 
5.8) the EM are overestimated in several pixels because 
of the same effects discussed in appendix B, because of 
the lack of constraints at lower temperatures. The noisy 
results at the low and high temperature end are expected 
considering the poor constraints provided by the selected 
EIS lines at those temperatures. 

With AIA, though some structures such as the side 
view loops are still apparent in the reconstructed EM 
maps to a lesser extent, the EM maps derived from the 



analysis of the synthetic data are not reproducing the 
true maps very accurately, especially for the top view 
case. For instance for the top view case of snapshot C 
(14) the maps derived from AIA present significant EM, 
especially in the emerging flux region, in the bin tem- 
perature log T[i^] ~ 6.4 — 6.7, where the real EM is zero 
almost everywhere. 

In Figure 20 we show, for all the cases including the 
effect of the noise, the scatter plots of the derived vs. 
true EM, for the central four temperature bins, i.e., 
log T[i^] : 5.5-5.8, 5.8-6.1, 6.1-6.4, 6.4-6.7. These plots 
show that, when integrating in broad temperature ranges 
(AlogT = 0.3), the EM derived from EIS reproduce 
rather well the true EM, while the EM derived from AIA 
have significantly larger scatter for the whole range of 
EM values. We note that these bins are significantly 
wider than typically used by observers when deriving 
EMDs. 

3.2. Full emission measure distributions 

In the previous section we have addressed a comparison 
of the general properties of the EMDs derived from syn- 
thetic EIS and AIA data with the input emission measure 
distributions from 3D rMHD simulations. In this section 
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Fig. 9. — Maps showing the Xq for the reconstruction of the EMD using the AIA intensities, for snapshot C. As in Figure 5, we show the 
results for the two LOS, "xy" (top), and "xz" (bottom), and also the case including the Poisson noise (right column). 



we will assess how well the derived EMDs reproduce the 
true EMDs in their details, on the finest temperature 
scale, and explore to what level of detail the MCMC 
EMD reconstruction method is reliable. 

In Figure 29 we have already shown the full EMDs for 
a selection of pixels chosen to explain the systematic dis- 
crepancies in T{EMDmax)- That set of plots already 
shows that, while as discussed above some general fea- 
tures of the EMD are recovered by the MCMC method, 
when comparing the EMDs on the fine temperature scale 
(i.e., here we used A log T[i^r] = 0.05 which is the typical 
intrinsic resolution of atomic data, and also the typical 
resolution used by observers), the derived EMD can sig- 
nificantly depart from the true distributions. 

In Figure 21 we show the comparison of the true and 
derived EMD for a pixel selected for each case (for every 
snapshot and LOS, for the cases including noise), to show 
examples where, even if the intensities are reproduced 
within the errors (x^ < 1), the derived and true EMD 
present significant discrepancies. 

First we note that in most of the examples of Figures 29 
and 21, there are several temperature bins within the 
logT = [5.5 — 6.7] range where the true EMD is not 
compatible with the derived EMD within the uncertain- 
ties. In general we find that the EMDs resulting from 
the analysis of EIS synthetic data reproduce the general 
properties of the EMDs - such as bulk of the emission 



and general shape - significantly better than the cor- 
responding AIA EMDs, as also shown in the previous 
section 3.1. However, even in the temperature range 
log T[i^] ~ [5.5 — 6.7] where the emission measure distri- 
bution should be well constrained by the selected spec- 
tral lines, the best fit EMD can be very noisy and present 
peaks and valleys which are not present in the true EMD. 

We note that, as discussed by Kashyap & Drake (1998) 
and Testa et al. (2011), the uncertainties associated with 
the EMD solution are correlated in the different temper- 
ature bins, and therefore, what are shown as error bars 
in the single temperature bins cannot strictly be inter- 
preted as error bars for the EMD value in that bin, and 
instead they describe the range containing 68% of the 
sets of solutions. Nevertheless, since they are typically 
used as error bars we investigate here also the limitations 
of this assumption. 

In order to assess the ability of the MCMC method 
to recover the full EMD for all the different cases stud- 
ied, we use two different ways to parametrize the "good- 
ness of the fit". First, we use the full information of the 
set of EMD solution produced by the MCMC routine 
for each pixel. For a given pixel, and a given temper- 
ature bin (we use here the full temperature resolution, 
i.e., A log T[i^] = 0.05), we can calculate the fraction of 
solutions included in the range between the true EMD 
value and the best fit value. This can be interpreted as 
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Fig. 10. — Xq maps analogous to Figure 9, for snapshot H. 

a probability for the EMD value to be part of the distri- 
bution of Monte Carlo solutions. Since it is not trivial 
to combine the results of the different temperature bins, 
we compute an average of these fraction over all temper- 
ature bins between logT = 5.5 and 6.7. If the average 
is close to 1, this implies that for most temperature bins 
the true EMD is far from the bulk of the distribution. 
We show the maps of the values obtained from the EIS 
EMDs in this fashion in the left panels of Figure 22 (for 
snapshot C) and 23 (for snapshot H). The corresponding 
maps for the AIA EMDs are not shown but are qualita- 
tively very similar to the EIS cases shown here. For the 
top view cases the maps present similar structuring to 
the Xo iiiaps shown in Figures 5 and 6, with the worst 
match of EMD found in the regions with large superpo- 
sitions of plasma structures with different densities and 
temperatures. For the side view, while lower in the atmo- 
sphere the values are similar to the top view case, they 
become ~ 1 everywhere at larger heights. The reason 
for this is that at larger heights the true EMD are close 
to isothermal, dropping to zero everywhere outside the 
peak region. For all the temperature bins where the true 
EMD is zero, the derived EMDs are non-zero causing 



Xo int 

0.0 0.5 1.0 1.5 2.0 




pix [0.6"] 



0.0 0.5 f.S' 1.5 2.0 




pIx [0.6"] 



large discrepancies in most of the temperature bins. 

Another way to calculate a measure of how well 
the derived EMDs match the true EMDs is to calcu- 
late a Xq{EMD) for the EMD, which we define as = 
Y.j[{EMDj^MCMC-EMDj^true)/(j{EMD)jf/df,^hem 
EMDj, 

true and EM Dj^MCMC ai"G respectively the val- 
ues of the true and derived EMD for the temperature bin 
j, (j{EMD)j is the uncertainty of the EMD in the T bin 
j defined by the range including 68% of the solutions, 
as described above, and df are the degrees of freedom 
(as defined at the beginning of section 3). We find that 
when calculating the so defined Xo{EMD) using the fine 
temperature grid (AlogT = 0.05) the Xo{EMD) values 
are very large {':^ 10) for both the EIS and AIA cases. 
Therefore we decrease the temperature resolution until 
the Xo{EMD) values become on average more reason- 
able. We find that this threshold is AlogT ~ 0.2 for 
EIS, and we show the corresponding maps of Xoi^MD) 
in the middle panels of Figure 22 and 23. These maps 
show the usual pattern of the features with large super- 
position of different structures, where the agreement is 
worse, but for the other regions the integration over large 
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Fig. 11. — Dependence of AIA temperature responses on plasma pressure (ue x T [cm~^ K]). In the pressure-temperature space, we 
plot maps of the ratio of the AIA response to the response calculated for a pressure of 10-*^^ cm""^ K, which is the pressure assumed for 
the standard AIA temperature responses available in SolarSoft. From left to right: top - 94A, 13lA, 17lA; bottom - 193A, 21lA, 335A 
passbands. We overplot the default AIA temperature responses (white lines; in units of [DN pix"-*^ cm^], and normalized to their 
peak), also showing the contribution of the lines (dashed lines) and the continuum emission (dot-dashed lines). 



temperature bins yields a better match of the derived 
and true EMDs, by smoothing out the large variabil- 
ity on small temperature scale shown in Figures 29 and 
21. We note that for the side view cases, high in the 
corona, the uncertainties of these EMD on the coarser 
temperature grid appear even overestimated, since the 
Xo{EMD) drops to very low values. This is consistent 
with the findings of Landi et al. (2012) who explored the 
ability of the MCMC methods to diagnose isothermal 
EMD. Landi et al. (2012) find indeed that for isother- 
mal plasma a bin width of AlogT ~ 0.05 is sufficient 
to diagnose the temperature distribution from spectral 
data. For AIA, even with AlogT - 0.3 the xl{EMD) 
values (maps are shown in the right panels of Figure 22 
and 23) are large for the top view cases, in agreement 
with the results discussed in above section 3.1. The side 
view cases show that with AIA the temperature binsize 
AlogT ~ 0.3 is adequate to diagnose isothermal plas- 
mas, as shown by the low Xoi^^^) values at high z 
values, but is still insufficient to guarantee a good match 
of the true EMD where the distributions are significantly 
multi-thermal. The EMDs derived AIA are a much less 
accurate reproduction of the true EMDs compared to 
EIS. This points to a significant problem in relying on 
AIA data exclusively, to diagnose the plasma tempera- 
ture distribution. This can be explained by considering 
two main factors: (1) the smaller number of constraints 
- AIA data provide a considerably more limited number 
of constraints (at best 6) compared to EIS; (2) imaging 



data have significantly more limited temperature diag- 
nostics than spectral data, because of their broader tem- 
perature sensitivity (even for the narrow AIA passbands) 
compared to the resolved spectral lines. In this respect, 
though outside the scope of this paper, we note that nar- 
row band coronal imagers also have the additional disad- 
vantage with respect to broad band instruments of being 
critically sensitive to the uncertainties in the atomic data 
(see, e.g.. Testa et al. 2012). We note that the EMDs 
derived from AIA present significant discrepancies with 
the true EMDs even in the case that does not include 
noise, or in the test case with very high S/N discussed at 
the beginning of this section when we discussed the fit 
to the observed intensities. We therefore conclude that 
the limitations in diagnosing EMDs with AIA data are 
mainly due to the intrinsic characteristics of the (broad) 
AIA temperature responses, and that the noise is not the 
main cause of the discrepancies between true and derived 
EMDs. 

While our results imply that in general the MCMC 
method does not reproduce the true EMDs at the highest 
temperature resolution, it is preferable for several reasons 
to run the method at high resolution (AlogT ^ 0.05) 
and then rebin the solution a posteriori. As discussed 
above (and shown in Figure 22 and 23), the high res- 
olution allows one to find nearly isothermal solutions. 
Also, we have run for a few pixels the MCMC on the 
EIS synthetic data using different values of A log T and 
found that adopting a too large bin size leads to less 
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Fig. 12. — Emission measure maps for the top view of snap- 
shot C, from the model (left panels) and from EIS synthetic data 
(without noise or including noise, in the middle and right panels 
respectively), integrated in six temperature bins: logT[X] : 5.0- 
5.5, 5.5-5.8, 5.8-6.1, 6.1-6.4, 6.4-6.7, 6.7-7. White corresponds to 
values of EM equal to zero. 

stable solutions, due to the rebinning of the temperature 
response function and consequent loss of temperature in- 
formation. For instance if the temperature bins become 
significantly large compared to the peaks in the tempera- 
ture response functions the rebinned responses often end 
up with peaks which are much less prominent, and sig- 
nificantly displaced in temperature. 

4. DISCUSSION AND CONCLUSIONS 

We have presented the results of a test of the limi- 
tations of plasma temperature diagnostics that are cur- 
rently available from spectral and imaging observations 
of the solar corona. Determining the plasma emission 
measure distribution is fraught with difficulties, both 
because of the intrinsic nature of the mathematical in- 
version problem, which is ill-posed, and because of the 
limitations of the available data. While imaging instru- 
ments provide only moderate temperature diagnostics 
(compared with spectrographs that can resolve spectral 
fines), they typically provide significantly better spatial 
coverage and resolution, and temporal cadence. 

Here we used advanced 3D radiative MHD simulations 
of the solar atmosphere (Hansteen et al. 2007; Gudiksen 
et al. 2011) as reahstic test cases, from which we pro- 
duced synthetic observables, and applied a Monte Carlo 
Markov chain forward modeling technique to derive the 
emission measure distributions (EMD). We then com- 
pared the results of EMD reconstruction from imaging 
(AIA) and spectral (EIS) synthetic data (based on the 
coronal properties in the models) with the input "true" 
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Fig. 13. — Comparison of emission measure maps obtained from 
the analysis of EIS synthetic data of snapshot C, analogous to 
Figure 12 but for the side view. 



EMD to establish the limitations of the temperature di- 
agnostic power of the available instruments, and how 
these results depend on the characteristics of the un- 
derlying thermal distributions. We also investigated the 
effect of the photon counting noise, by running the analy- 
sis with or without randomization accounting for Poisson 
noise. 

These 3D simulations provide us with the opportunity 
to improve upon previous work by exploring more real- 
istic configurations, with significant superposition of dif- 
ferent structures along the LOS, allowing a statistical ap- 
proach to determine the accuracy and limitations of the 
plasma diagnostics, for a variety of spatial and thermal 
structuring of the plasma. The three dimensional nature 
of the simulations also allows us to explore a variety of 
realistic viewing angles, reproducing typical distributions 
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Fig. 14. — Comparison of emission measure maps obtained from 
the analysis of the top view AIA synthetic data of snapshot C, 
analogous to Figure 12. 



of structuring ranging from on disk to limb observations. 

We assess the robustness of the EMD reconstruction 
method by using several parameters. We explored the 
ability of the method in: (a) reproducing the "measured" 
intensities; (b) determining the temperature of the peak 
of the emission measure distribution; (c) deriving the 
emission measure values in broad temperature bins; (d) 
reproducing the true emission distribution in its details. 

The analysis of the spectral EIS synthetic data show 
that the measured intensities are reproduced generally 
well by the inferred EMD, even when including the ef- 
fect of noise, with the exception of regions with mixing 
of regions with significantly different densities along the 
LOS. The temperature of the peak of the EMD is repro- 
duced reasonably well for the side view, in larger areas 
where the EMD is close to isothermal. For the top view, 
where the EMD are broad and in particular where there 
are large amounts of dense material and large superposi- 
tion of different structures (emerging fiux region, moss), 
the temperature of the peak EMD appears to be system- 
atically underestimated. The maps of emission measure, 
when integrated in broad temperature bins, are similar 
to the true distributions, in the well constrained range 
(log(T[i^]) ~ [5.6 — 6.4]); the low and high temperature 
ends of the considered range are rather noisy and present 
spurious components. However, when considering the 
direct comparison of the derived EMD curves with the 
true distributions, at the full temperature resolution, we 
find that the detailed properties of the true EMD on 
the smallest temperature scale are not accurately repro- 
duced. The best fit EMD can present peaks which are 
not displayed by the true EMD, or vice versa, narrow 




Fig. 15. — Comparison of emission measure maps obtained from 
the analysis of the side view AIA synthetic data of snapshot C, 
analogous to Figure 12. 



peaks might be missed by the reconstruction method. 
The worst cases show how the uncertainties do not gen- 
erally account for these discrepancies, and are therefore 
likely underestimated by the method. We explore the 
temperature resolution at which the discrepancies be- 
tween true and derived EMD become smaller, and we 
find that for EIS the AlogT needs to be > 0.2. We 
therefore conclude that the temperature diagnostic pro- 
vided by spectral data with enough observational con- 
strains and good temperature coverage, are on average 
reliable to derive the general characteristics of the emis- 
sion measure distribution, such as the emission measure 
in broad temperature bins, but that the results are not 
robust enough as far as the fine details of the EMD, such 
as the presence of narrow peaks, are concerned, and that 
features on scales below > 0.2 in logT cannot be trusted. 
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Fig. 16. — Emission measure maps for the top view of snap- 
shot H, from the model (left panels) and from EIS synthetic data 
(without noise or including noise, in the middle and right panels 
respectively), integrated in six temperature bins: \ogT[K] : 5.0- 
5.5, 5.5-5.8, 5.8-6.1, 6.1-6.4, 6.4-6.7, 6.7-7. Figure analogous to 
Figure 12, but for snapshot H. 

The results obtained from the analysis of the AIA syn- 
thetic data are more unsatisfactory. The synthetic in- 
tensities for the side view LOS are low in some chan- 
nels, making the results more sensitive to the effect of 
the noise. The maps of temperature at which the EMDs 
peak show that, as for the EIS case, the temperatures are 
not very well reproduced, and appear systematically un- 
derestimated, especially in the top view case. Also the 
maps of emission measure integrated in broad temper- 
ature bins are a poor match for the true distributions, 
especially for the top view case. Finally, the direct com- 
parison of derived and true EMD curves at the full tem- 
perature resolution indicates discrepancies well beyond 
the errors, implying that the EMD are not at all well con- 
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Fig. 17. — Comparison of emission measure maps obtained from 
the analysis of the side view EIS synthetic data of snapshot H, 
analogous to Figure 16. 

strained by the AIA data, often not even in their general 
characteristics. Even degrading the temperature resolu- 
tion to AlogT = 0.3 the derived EMD do not match 
accurately the true EMD, especially where the tempera- 
ture distributions are broad. We therefore conclude that 
the limited number of constraints and the broad tem- 
perature sensitivity of the AIA passbands, which often 
present multiple peaks, critically hamper the capabil- 
ity of diagnosing the plasma temperature distribution 
on the basis of AIA imaging observations exclusively. 
The addition of simultaneous Hinode/XRT data, which 
provide complementary imaging observations in X-ray 
broadbands, might improve considerably the ability to 
reconstruct EMDs from AIA data, by providing addi- 
tional constraints, especially at the high temperature end 
(logT ^ 6.5 — 7), where AIA provides limited informa- 
tion. We plan to explore this issue in a follow-up paper. 
In this paper we have not addressed the use of EIS and 
AIA together because, as discussed here above and in 
other papers (see e.g.. Warren et al. 2011; Testa et al. 
2011), the spectral data, when available, provide much 
better constraints to the temperature distribution of the 
plasma, compared with imaging data, and the imaging 
data provide only limited additional constraints. How- 
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Fig. 18. — Comparison of emission measure maps obtained from 
the analysis of the top view AIA synthetic data of snapshot H, 
analogous to Figure 16. 



ever, imaging data provide a much better temporal and 
spatial coverage of the coronal plasma, so we focused on 
investigating the usefulness of imaging data for thermal 
diagnostics, when spectral data are not available, as in 
the large majority of solar observations. 

The approach adopted here, using 3D simulations, has 
provided us with the opportunity to improve upon previ- 
ous work by exploring more realistic configurations, with 
significant superposition of different structures along the 
LOS, allowing a statistical approach to determine the 
accuracy and hmitations of the plasma diagnostics, for a 
variety of spatial and thermal structuring of the plasma. 
These results provide a stringent and accurate assess- 
ment of the limitations of the temperature diagnostics, 
indicating the extent of the sensitivity of EMD measure- 
ments to the structuring along the LOS, and to char- 




log, J:,6,70-7.00 tri 



Fig. 19. — Comparison of emission measure maps obtained from 
the analysis of the side view AIA synthetic data of snapshot H, 
analogous to Figure 16. 

acteristics inherent to the instrumentation. However, we 
emphasize that the exercise we carried out focuses on the 
use of the MCMC reconstruction method, and therefore 
it only assesses the robustness of this particular method, 
which is however largely adopted and trusted for stud- 
ies of emission measure distributions of solar and stel- 
lar coronal plasmas. We note that these results actually 
present an optimistic view, because they do not account 
for several effects that further complicate real observa- 
tions, such as for instance uncertainties in both instru- 
ment calibration and atomic data, infiuence of blends 
in fitting the spectral data, and distribution of element 
abundances. Finally, the 3D simulations adopted here 
model limited coronal volumes, therefore likely underes- 
timating the effects of superposition along the LOS. 
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Fig. 20. — Scatter plots of EM derived from EIS (first two columns) and AIA (right columns) synthetic data from the two snapshots 
(cool simulation "C" in the top row, and hot simulation "H" in the bottom row), vs. the corresponding true EM values. For each of the 
panels, we plot the values for the central four temperature bins (\ogT[K] : [5. 5-5. 8], [5. 8-6.1], [6. 1-6. 4], [6. 4-6. 7]), and excluded the EM values 
for the two temperature bins at the extremes of the range (\ogT[K] : [5. 0-5. 5], [6. 7-7.0]). 



:(x,y)=13,^ 
10^^ rXo'int=0.75 




:(x,y) = 20,1 1 
10^^ rXo'int=0.76 





:(x,z) = 4,14 
rXo'int=0.70 





5.0 5.5 6.0 6.5 7S).0 5.5 6.0 6.5 7S).0 5.5 6.0 6.5 7S).0 5.5 6.0 6.5 7.0 
Log T[K] Log T[K] Log T[K] Log T[K] 



Fig. 21. — Sample of comparisons of true (red curves) and derived EMD (including noise), for one pixel for each case. The results from 
EIS synthetic data are shown in the left two columns, and the ones from AIA in the right two columns. The results for snapshot C are 
shown in the top row, while the bottom four panels show the results for snapshot H. 
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Fig. 22. — Left: Maps of the "goodness of fit" for the EMD for the EIS case both top view (top row) and side view (bottom row) for 
snapshot C. The value in each pixel represents the average over all temperature bins between logT = 5.5 and 6.7, of the fraction of the 
solutions falling in the range between the best fit EMD and the true EMD which is used as a measure of the probability that the true EMD 
is part of the same distribution of the EMD solutions. A value close to 1 implies that for most temperature bins the true EMD is far from the 
bulk of the distribution. We use the full temperature resolution (AlogT[X] = 0.05). Middle: Maps of x^{EMD) (see text for definition), 
parametrizing the ability of the MCMC method to recover the full EMD in each temperature bin from EIS data (including noise). We used 
Alogr[K] = 0.2. Right: Maps of x^{EMD) for the EMD derived from AIA data (including noise), and assuming AlogT[X] = 0.3. 



were discussed with other colleagues. 
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Fig. 23. — Plots analogous to plots in Figure 22 but for snapshot H. 
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Fig. 24. — SDO/AIA 193A channel full disk observation taken on 2012 January 6 around 23UT, where we selected a subset of small 
regions (100 x 100 pixels) of different coronal features: quiet Sun (QS), cool fan loops (CL), limb (LM), and active region plasma (ARl, 
AR2). In Figure 25 we show the distributions of the 6 AIA EUV coronal channels intensities for these selected regions, to be compared 
with the distributions of the synthetic intensities derived from the simulations (shown in Figure 26). 



APPENDIX 

A. COMPARISON OF SYNTHETIC INTENSITIES WITH ACTUAL OBSERVATIONS 

In this appendix we present a comparison of the synthetic intensities derived from the simulations, as described in 
section 2.1, with measured intensities from recent SDO/AIA observations. We considered AIA observations taken on 
2012 January 6, around 23UT (the 193 A full disk image is shown in Figure 24), and selected a few small regions (100 
pixels X 100 pixels, corresponding to ~ 60 arcsec x 60 arcsec) sampling a variety of coronal features, ranging from 
quiet Sun to active regions, and including areas both on-disk and above limb. 

In Figure 25 we show the distribution of the observed intensities (in units of DN s~^ pix~^) in the six selected AIA 
EUV bands, for the different regions shown and labeled in Figure 24. These plots clearly show that: (a) average 
absolute intensities change by more than an order of magnitude from region to region; (b) the distribution in each 
channel, in each of the regions is rather narrow; (c) the distributions of intensities in the 94A, 131A, and 335A channels 
typically peak at values about two orders of magnitude lower than the intensities in the other three bright channels; 
(d) the intensities in the 211 A, and 335 A channels typically have relative lower values in cooler regions (quiet Sun and 
cool fan loops; regions QS and CL). 

In Figure 26 we show the distribution of the synthetic intensities in the AIA passbands, for the two lines of sight for 
each of the two snapshots. The smaller, cooler snapshot (C) is characterized by intensities similar to the values observed 
in cool fan loops (region CL) and quiet Sun (region QS), whereas the larger and hotter snapshot (H) produces intensities 
more similar to the observed values of active region plasma (regions ARl and AR2). We note that, as expected, the 
top view ("xy") provides a more realistic comparison for real coronal observations, while for the side view ("xz") the 
distributions have tails at low intensity values, which are not observed in the AIA selected data. This is because the 
small box size does not accurately capture the enormous line-of-sight at the limb in the optically thin corona. However, 
for the side view of snapshot C the inversion of the peaks of the 335A and 131A channels with respect to the top 
view (i.e., for the side view the 131 A distribution peaks at lower intensity values than the 335 A) reproduces what is 
observed in actual data at the limb (region LM) compared with quiet on disk regions (QS, CL). We also note that the 
distributions for snapshot H appear generally broader than the observed distributions. For snapshot C the distribution 
of each passband peaks around a different value, whereas the observed distributions typically show a common peak 
for the weaker channels and a separate common peak for the stronger channels. Finally, the intensities in the 94 A 
channel are systematically lower than the observed values, but this is expected on the basis of recent studies arguing 
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Fig. 25. — Histograms showing the distributions of the observed intensities (in DN s"-*^ pix"-*^) in the AIA channels (shown in different 
colors, as labeled in the top left panel, QS), for the regions indicated in Figure 24: quiet Sun (QS), cool fan loops (CL), limb (LM), and 
active region plasma (ARl, AR2). 

for incompleteness of atomic databases, especially in the 94A passband (e.g., Testa et al. 2012). By and large, these 
qualitative comparisons indicate that the simulations produce emission values of similar order of magnitude of real 
observations, and therefore represent a reasonable test case for the coronal diagnostics. 

B. EMISSION MEASURE PEAK TEMPERATURE 
See complete preprint version at http://folk.uio.no/bdp/papers/3dEMD_ptesta.pdf 
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Fig. 26. — Histograms showing the distributions of the AIA intensities synthesized from the two snapshots (left panels: smaller cooler 
snapshot, C; right panels: larger hotter snapshot, H), and for the two diflFerent views ("xy" in the top panels and "xz" in the bottom 
panels). 



